Electron tunneling into a quantum wire in the Fabry-Perot regime 
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We study a gated quantum wire contacted to source and drain electrodes in the Fabry-Perot 
regime. The wire is also coupled to a third terminal (tip), and we allow for an asymmetry of the tip 
tunneling amplitudes of right and left moving electrons. We analyze configurations where the tip 
acts as an electron injector or as a voltage-probe, and show that the transport properties of this 
three-terminal set-up exhibit very rich physical behavior. For a non-interacting wire we find that a 
tip in the voltage-probe configuration affects the source-drain transport in different ways, namely by 
suppressing the conductance, by modulating the Fabry-Perot oscillations, and by reducing their vis- 
ibility. The combined effect of electron electron interaction and finite length of the wire, accounted 
for by the inhomogeneous Luttinger liquid model, leads to significantly modified predictions as com- 
pared to models based on infinite wires. We show that when the tip injects electrons asymmetrically 
the charge fractionalization induced by interaction cannot be inferred from the asymmetry of the 
currents flowing in source and drain. Nevertheless interaction effects are visible as oscillations in 
the non-linear tip-source and tip-drain conductances. Important differences with respect to a two- 
terminal set-up emerge, suggesting new strategies for the experimental investigation of Luttinger 
liquid behavior. 

PACS numbers: 73.23.-b, 71.10.Pm, 73. 23. Ad, 73.40.Gk 



I. INTRODUCTION 

Electron scanning of a conductor with a probe ter- 
minal is a customary technique to investigate its local 
properties. The local density of states can be gained 
from the dependence of the tunneling current on the ap- 
plied bias. Nowadays, atomically resolved images arc 
obtained both with scanning tunnel microscopes (STM) 
and atomic force microscopes (AFM)ii So far, most of 
the efforts of the scientific community have focused on 
improving the resolution power of the probe terminal. 
For instance, the recent realization of stable and sharp 
superconducting STM tips exploits the singularity in the 
quasiparticlc density of states to this purposci^ A probe 
terminal, however, may also be used as a "handle", i.e. 
as an active component to tune the transport properties 
of the conductor. Recent works in this direction have 
shown that the sign of the supercurrent can be changed 
when a third terminal injects electrons into a Josephson 
junction under appropriate conditions;^ that the conduc- 
tance of a quantum dot can be tuned by moving an AFM 
tip over the sample^ or that a single-electron transistor 
can be used to cool down a nanomechanical resonator, or 
to drive it into a squeezed statci^ 

The promising applications of scanning probes in the 
study of transport properties of nanodevices require a 
theoretical analysis of electron transport in a three- 
terminal set-up, a subject which has been explored only 
partly so far. In particular, most of the available in- 
vestigations are restricted to the case of non-interacting 
conductors whereas relatively little attention has been 
devoted to those nanodevices in which electronic corre- 
lations play a dominant role. This is the case for one- 



dimensional (ID) conductors, such as semiconductor het- 
erostructurc quantum wires^ and single-walled carbon 
nanotubesi^i^ There, electron-electron interaction dra- 
matically affects the dynamics of charge injection. The 
response of the system to the scanning probe is quite 
different from that of ordinary three dimensional met- 
als, since in ID electronic correlations lead to a break- 
down of the Fermi liquid picture. Semiconductor quan- 
tum wires and carbon nanotubes rather exhibit Luttinger 
liquid (LL) behavior j ^^i^^i^^i^"^ While for this type of sys- 
tems two-terminal electron transport has been widely an- 
alyzed in the last 15 yearS j^'^i-'^'-'i^-'^i-'^^'-'^'^i"'^^ the electric cur- 
rent and noise in a three-terminal set-up, including source 
and drain electrodes and a tip, have remained mostly un- 
explored. 

There are, however, a few notable works in this di- 
rection. The case where a bias is applied between a tip 
and a semi- infinite LL was investigated by Eggert^ and 
by Ussishkin and Glazmaufi^ Martin and co-workersi2ii^ 
have recently analyzed the electric noise of the current in- 
jected from a tip into a nanotube adiabatically contacted 
at each end to grounded metallic leads. 

In this paper we extend these investigations to a quite 
general three-terminal set-up. We shall thus explore the 
non-equilibrium current in all three terminals in pres- 
ence of a transport voltage between the source and drain 
electrodes, an applied tip voltage, and also a tunable 
gate voltage. This enables us to address various physical 
phenomena that are of relevance for recent experiments. 
Among other effects, we discuss the influence of the tip 
on the transport along the interacting wire, even when 
no net current is injected from the tip into the wire. In 
particular, we focus on the Fabry-Perot transport regime 
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FIG. 1: (Color online) Sketch of the set-up. A quantum wire 
is connected to two metallic electrodes, denoted as source (S) 
and drain (D) at voltages Vs and Vd, respectively. A third 
sharp electrode, denoted as tip (T), at voltage Vt, injects elec- 
trons into the wire at position xq. A gate (G) is also present 
and held at a gate bias voltage Vg- The contact resistances 
are accounted for by two delta-like scatterers with strengths 
Ai and A2 , and the electron tunneling amplitudes between the 
tip and the wire are denoted by 7±. 



of the wire, which could be recently observed in carbon 
nanotubeaiS'i^ii^'SS, and analyze how Fabry-Perot oscil- 
lations are modified by both the presence of the tip and 
the electron-electron interaction. To this purpose, the 
finite length of the wire, the contact resistances at the 
interfaces between the wire and the side electrodes, as 
well as an arbitrary position of the tip along the ID wire 
are taken into account in our model. Furthermore, in- 
spired by recent experiments on semiconductor quantum 
wireai^"^, we allow for an asymmetry in electron tunnel- 
ing from the tip, and investigate how the presence of side 
electrodes affects the fractionalization of charges injected 
by the tip into an interacting wire. Finally, regarding 
the experimental observation of interaction effects, we 
discuss the advantages of a three-terminal set-up over a 
two-terminal one. 

The paper is organized as follows. In Sec. |TT] we de- 
scribe the model that we adopt for the set-up. In Sec. IIIII 
we provide results about the electric current in the case 
of a non-interacting wire, while Sec. llVl is devoted to the 
effects of electron-electron interaction. Finally, we shall 
discuss the results in Sec. IVland present our conclusions. 
Some more technical details are given in the appendices. 



II. THE MODEL 



the injection point on the tip. The latter is located 
at position xq with respect to the middle of the wire, 
and electron injection is modeled by a tunnel Hamilto- 
nian. We also envisage the presence of a metallic gate 
(G), biased at a voltage Vq. Screening by this gate 
yields an effectively short-ranged electron-electron inter- 
action potential within the wire, for which the LL model 
applies i^i^i^i The total Hamiltonian of the system reads 



(1) 



where the first term describes the wire and its coupling to 
the S and D electrodes as well as to the gate. The second 
term accounts for the tip, and the last one describes wire- 
tip tunneling. 

As far as the wire is concerned, we shall address here 
the low-energy regime, where the wire electron band can 
be linearized around the Fermi level. Then the wire elec- 
tron operator '^(x) can be decomposed into right- and 
left-moving components '^-^-(x) and ^'_(a;) 
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where few denotes the equilibrium Fermi momentum of 
the wire. By definition, this is the Fermi momentum 
in case that the electrochemical potentials of all elec- 
trodes, source, drain, tip and gate, are identical. This 
corresponds to vanishing applied voltages. Explicitly the 
Hamiltonian of the wire reads 



In Eq. ^ the first term 
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describes the band energy linearized around the wire 
Fermi points ±fcw and characterized by a Fermi velocity 
wyv- The symbol : : stands for normal ordering with re- 
spect to the equilibrium ground state. The second term 
models scatterers at the interfaces^i^ with the S and D 
electrodes 



(5) 



We consider a single channel spinless quantum wire 
connected, as sketched in Fig. [Tl to two metallic elec- 
trodes, source (S) and drain (D), as well as to a third 
sharp electrode, henceforth denoted as tip (T). The wire 
has a finite length L and for the x coordinate along it 
we choose the origin in the middle of the wire so that 
the interfaces to the S and D electrodes are located at 
xi = —L/2 and X2 = +Ll'i^ respectively. Electron 
backscattering at the side contacts due to non-adiabatic 
coupling is modeled by two delta-like scatterers. The tip 
is described as a semi-infinite non-interacting Fermi liq- 
uid, and y < denotes the coordinate axis along the 
tip orthogonal to the wire, the origin corresponding to 



where the dimensionless parameters > denote the 
impurity strengths at the contacts x^, and the term p{x)=: 
■.'^\x)'^ {x): is the electron density fluctuation with re- 
spect to the equilibrium value. The third term in Eq. ([3]), 
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with 



pxv {x) p{x) dx (6) 



for x < -L/2 

for -L/2<x< L/2 , (7) 
for X > L/2 
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accounts for the bias Vs and Vd of the source and drain 
electrodes, as well as for the gate voltage Vg- The applied 
transport voltage is then V — Vs --Vb- Finally, the last 
term 
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/ dx V 



: pr{x)pr'{x) 



(8) 



describes the screened Coulomb interaction in the 
wire^^i^^ where pr{x)= :^J(a;)5'^(a;): is the density fluc- 
tuation of r-moving electrons. As it is customary in 
LL theory, in the sequel, we characterize the interaction 
strength by the dimcnsionlcss coupling constant 



(9) 



The Hamiltonian of the tip, the second term in Eq. (fT|). 
reads 



Here 



Ht — Hkin.T + Hn 



ihvT I dy : c\y)dyc{y) : (11) 



(10) 



describes the (linearized) band energy with respect to the 
equilibrium Fermi points ifcx of the tip, and vt denotes 
the Fermi velocity. Notice that the integral runs also over 
the positive y-axis, since right and left moving electron 
operators along the physical tip axis y < have been 
unfolded into one chiral (right-moving) operator c{y) de- 
fined on the whole y-axis. The second term in Eq. (fTO|) 
describes the bias Vt applied to the tip which affects the 
incoming electrons according to 



eVr / dy : c\y)c{y) 
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Finally, the third term in Eq. Ilj accounts for the wire-tip 
electron tunneling and reads 



Htun = nV^^ J2 > (e''"'^""°*];(xo)c(0) + h.c.) 



(13) 

where 7^ is the dimensionless tunneling amplitude for r- 
moving electrons, and x = xq [y = 0] is the coordinate 
of the injection point along the wire [tip]. Here, we have 
allowed for a right/left asymmetry of electron tunneling 
between the tip and the wire, which can arise from the 
presence of a magnetic field^^'^'^l. Note that for 7+ 7^ 
7_ the Hamiltonian is not invariant under time-reversal 
symmetry. 

In the following sections the electron current will be 
evaluated in the three terminals of the described set-up. 
Explicitly we shall compute 

I{x,t) = evw (14) 
(: '^l{x,t)^+ix,t) : - : *L(x,t)«'_(x,t) :) 



where x, with > L/2, is a measurement point located 
in the S or D leads. As far as the tip is concerned, due 
to the unfolding procedure described above, the electron 
current flowing in the tip at a point y < acquires the 
form 

I{y,t) = evT (: {y,t)c{y,t) : - : c\-y,t)c{-y,t) :) . 

(15) 

In Eqs. (fH|) and the averages arc computed with 
respect to the stationary state in presence of the applied 
dc voltages Vs, Vd, Vt and Vq. 

Under these conditions, the current in each electrode is 
actually independent of the measurement point. We thus 
denote by /s and /d the currents flowing in the source and 
drain electrodes. The current /g is positive when flowing 
into the wire, while /d is positive when flowing out of the 
wire. The current /t flowing in the tip is positive when 
flowing in the direction of the tip- wire tunnel contact. 
Current conservation then implies Id = Is + It, so that 
all currents can be expressed in terms of two independent 
quantities. One can write 



h 
Id 



Im 
Im 



h/2 
It/2 



(16) 
(17) 



where Im describes the current flowing in the wire under 
the condition that no net current flows through the tip 
(voltage probe conflguration) . Importantly, Im should 
not be identified with the two-terminal current flowing in 
the absence of the tip. Indeed, while 7± = implies that 
It = 0, the opposite does not hold, so that Im needs to 
be evaluated by accounting for the whole three-terminal 
set-up. 



III. THE NON-INTERACTING CASE 

In this section we flrst discuss results for the case that 
the electron interaction ([8]) is neglected. Then the Hamil- 
tonian ([1]) of the whole system is quadratic in the fields 
'i!±{x) and c{y), and transport properties can be deter- 
mined within the Landauer-Biittiker formalism. In the 
three-terminal set-up that we are considering, the scat- 
tering matrix S{E) is a 3 x 3 matrix which depends on 
the energy E measured with respect to the equilibrium 
wire-lead Fermi level. The currents Im and It defined 
through Eq. HH) and ^ read 



^ / {\S,2\' + \S2in{fs{E)-ME))dE 
1 1"°° 

^/ M {.fs{E) - fT{E)) dE (18) 

~^ J —OQ 
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(ME) - fu{E)) dE 
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and 



/oo 
dE [\S3,\' (ME) - fs{E)) 
-OO 

+ \S32\HME) ~ ME))] . 



(19) 



In the S-matrix elements appearing in Eqs. (jlSp and ([T 
the source, drain and tip electrodes are identified as 1 



J 



2, and 3 respectively, whereas their Fermi functions are 
denoted as /s, /d and /t- Note that the S-matrix is in 
general not symmetric, because time-reversal symmetry 
is broken for 7-|_ ^ 7_. The S-matrix can straightfor- 
wardly be evaluated with standard techniques by com- 
bining the transfer matrices M^. (i — 1,2) of the two 
side contacts 
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with the one, M^^, at the tip injection point 
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(20) 
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Here, we have introduced the ballistic frequency 

v-p 

OJL^ -r 



(22) 



associated with the length of the wire, and the following 
dimensionlcss quantities 



^° = T' 
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The scattering matrix is obtained as a combination of the 
elements 
the form 



elements of the transmission matrix M = Mj-gMj-gMj-j in 
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- M21M32 
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M33^ 


/I22 - M23M32/ 



where My are the matrix elements of M. 



A. Fabry-Perot oscillations in a two-terminal 
set-up 

Before discussing the influence of the STM tip, we 
shortly describe the transport properties in the absence 
of the tip, i.e. for 7± = 0. In this case we have a two- 
terminal set-up with /t = and Is = Id = Im- The 
solid line in Fig. [2] shows the two-terminal conductance 
dlM./dV at zero temperature plotted in units of /h as 
a function of the (dimensionless) source-drain bias 



eiVs-Vp) 

hujL 



(25) 



for identical contact impurity strengths Ai = A2. 
For Ai <C 1 the conductance shows the typical 
Fabry-Pcrot oscillations with maximum values close to 
one. For carbon nanotubcs the Fabry-Perot regime of 
highly transparent contacts could be reached experimen- 
tally only recently due to progress achieved in device 
contactingi ^°i^^i^^i20 In the sequel, we will focus on this 
regime. 

The electron current Is = Id = Im can be written as 



= Iq 



'imp 



(26) 



where Jo = {e'^/h)V represents the current of a perfectly 
contacted wire, and /imp characterizes the (negative) cor- 
rection due to the contact resistances. The exact expres- 
sion for Jimp, which can be gained from the S-matrix, is 
not easily tractable for arbitrary impurity strengths and 
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temperature. In the Fabry-Perot regime at zero tem- 
perature, however, a simpler expression is obtained by 
expanding in terms of the impurity strengths. To third 
order in the A,;'s one obtains 



/ir 



27r 



(jinc H~ Jcoh) 



(27) 



where jinc and jcoh are dimensionless quantities describ- 
ing the incoherent and coherent contributions, respec- 
tively, to the reduction of the current by the contact im- 
purities. The term 



Jin 



(28) 



is linear in the applied bias voltage, and the coefficient 
of proportionality is the "classical" series resistance of 
two impurities. In contrast, the term jcoh stems from 
quantum interference between scattering processes. This 
interference leads to the Fabry-Pcrot oscillations of jcoh- 
Explicitly, 



Jcoh 



. ,-(2) , J3) 
J coll ' J coh 



(29) 



where 



,•(2) 
./coh 



-2A1A2 cos [2(7iw + 'tw — uq)] sin(u) (30) 



and 

jfl -2 (AiA^ + \l\2) sin [2(uw + kw - ug)] sin(u) , 

(31) 

where we have introduced 



WW 



2hujr, 



(32) 



From Eqs. (pn)) - ((5T|) one can see that Fabry-Perot oscil- 
lations arise both as a function of the source-drain bias 
u and as a function of the gate voltage uq. Note that for 
a non-interacting system the period in the former case is 
twice as large as the period in the latter case. 

We also emphasize that j^^^ originates from impurity 
forward-scattering processes (more precisely from second 
order in backward scattering and first order in forward 
scattering). Forward scattering processes are typically 
neglected in single impurity problems, where they can be 
gauged away. However, when two or more impurities are 

present they affect the coherent part of transport. Al- 

(2) 

though this contribution is in general smaller than j^^f^, 
it becomes the dominant term for the Fabry-Perot oscil- 



lations when 



.(2) 
J coh 



vanishes, which is the case for 



ejVs + Vu- 2Vg) 

2hjjjr 



odd integer . (33) 



Thus, the third order term is crucial for certain values of 
the biasing voltage. 

We conclude the discussion of the two-terminal case 
by emphasizing that for a non-interacting wire in the 



Fabry-Perot regime the current depends not only on the 
difference Vs — Vd, but in general on Vs and Vd sepa- 
rately. This is simply due to the fact that Fabry-Perot 
interference effects lead to an energy-dependent trans- 
mission coefficient and, hence, to non-linearity in the ap- 
plied bias. Notice that Eqs. jM]), jSO]) and jST]) fulfill the 
gauge-invariancc condition emphasized by Biittiker;^ 
since they are invariant under an overall shift of the po- 
tentials Vp Vp + const [p = S, D, G). 



B. Effect of the tip on Fabry-Perot oscillations 

In this section we shall address, within the non- 
interacting electron approximation, the effect of the STM 
tip on the Fabry-Perot oscillations. When 7± 7^ 0, the 
currents Im and It are non-vanishing for arbitrary val- 
ues of the applied voltages Vs, Vd and Vr- We analyze 
the effects of the tip as a function of the total tunneling 
strength 7, defined through 



„2 7i + 7^ 



the tunneling asymmetry coefficient 



7+ -71 
7++7- 



Ixl < 1 , 



(34) 



(35) 



and the position .tq of the tip. 



We start by considering the situation where the tip be- 
haves as an electron injector: a bias is applied between 
the tip and the source and drain electrodes, which, for 
simplicity, are assumed to be at the same electrochem- 
ical potential. A quite standard calculation applies to 
the case of fully symmetric tunneling (x = 0), allowing, 
e.g., to relate the local density of states in the wire to 
the non-linear conductance as a function of the tip-wire 
bias. Here, we shall instead focus on the case of fully 
asymmetric tunneling (x = ±1), which has become of 
particular interest due to recent experiments where only 
right-moving and/or only left-moving electrons could be 
selectively tunneled into a semiconductor quantum wire 
due to the presence of a magnetic field normal to the 
plane of the wire and the tip.— We find that novel phys- 
ical aspects emerge from a tunneling asymmetry. In the 
first instance, a direct inspection of the scattering matrix 
(j24p shows that its elements Sy are independent of xq, 
implying that, differently from the case of symmetric tun- 
neling X = 0, the lead currents /d and la, do not depend 
on the position of the tip. Furthermore, asymmetric tun- 
neling can be used to extract the transmission coefficient 
of each contact. Indeed evaluating the asymmetry 



A{x) 



|/d|-|/sI 



(36) 



between Jd and /s in the two cases of totally asymmetric 
injection only to the right (x = 1) and only to the left 
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Symmetric tunneling 





FIG. 2: (Color online) Zero temperature differential con- 
ductance as a function of the source-drain bias for a non- 
interacting wire characterized by contact impurity strengths 
Ai = A2 = 0.1 and a Fermi wavevector kw ~ 0.3. The tip 
is located in the middle of the wire and the tip voltage Vt is 
adjusted to fulfill the condition It = 0. Tunneling is symmet- 
ric (x = 0) and the tunneling strength has the values: 7 = 
(solid line), 7 — 0.1 (dashed-line) , 7 = 0.5 (dashed-dotted 
line) and 7 = 1 (dotted line). The gate is grounded (Vg ~ 0), 
and the bias is applied symmetrically (Vs/d = i''^/2). 



(x — ^1), one obtains 



and 



A+ = 



A- = 



1 + A?-Ai 



A? 



1 + Al 



Ai 



A? 



1 + xi 



Al 



(37) 



(38) 



where A± = ztA{±l). From these coefficients it is 
straightforward to extract the strengths of the contact 
impurities 



A? 



Xi = 



A+ +A^ 



1-A+ 
A+ + A- 



(39) 



as well as the transmission coefficients 

T ■ 1 _ A+ + A- 



(40) 



related to each of the two contacts. 

Notice that, while /g and /d depend on the tempera- 
ture T, Eqs. ([571) ^iid ([55)1 are independent of T within 



the approximation of a linearized band. Interestingly, 
these equations also enable one to identify the relation 
between the current asymmetry coefficients A± and the 
two-terminal conductance G2t = 9/m/3F|^=o- In Ref. 
[27l . the equality A± = G'2t/(e^//i) is claimed to hold for 
a set-up with symmetric contacts to the leads, even in 
the presence of interactions. However, Eqs. ([37]) and ((38l) 
show that for a quantum wire in the Fabry-Pcrot regime, 
even in the absence of interactions and with perfectly 
symmetric contacts Ai = A2, one has 



A+=A-^ G2t/{eVh) 



(41) 



since A± — 1/(1 -f- 2Xf) is a constant, whereas G2t 
depends on temperature, source-drain bias and gate 
voltage. The equality sign in Eq. (|¥T|l holds only under 
the specific circumstances of perfectly transmitting 
contacts (Ai,2 = 0), or of a perfectly symmetric set- 
up (Ai = A2 7^ 0) at sufficiently high temperatures 
fcer 3> fi^L, where Fabry-Pcrot oscillations of G2t wash 
out. 

The second situation that we want to investigate is 
when the tip voltage Vr is set to an appropriate value 
Vt so that no net current flows through the tip. This 
corresponds to a situation where the tip behaves as a 
voltage probe^^. Notice that, even under the condition 
It = 0, electrons can tunnel from the tip to the wire and 
vice versa, and therefore the tip does affect the electron 
transport between source and drain. 

We start by describing the case of symmetric tunnel- 
ing (x = 0) with the tip located in the middle of the wire 
{xq = 0). The differential conductance dlM_/dV, evalu- 
ated under the condition /x = 0, is depicted in Fig. [2] as a 
function of the source-drain bias ((25l) . for different values 
of 7, ranging from weak to strong tunneling. The tip has 
three main effects on the Fabry-Perot oscillations: i) an 
overall suppression of the conductance, ii) a modulation 
of the maxima and minima, and iii) a reduction of the 
visibility of the oscillations. 

The origin of the first effect can be illustrated already 
in the case of a clean wire (A^ = 0), where it is easy 
to show that the condition /t = is fulfilled for a tip 
voltage Vt = (Vs + Vd)/2, and that 



e ' Vs - Vu 
h I-H7V2 



^d), 



(42) 



Notice that a reduction of the conductance already shows 
up to order 7^ in the tunneling strength. The reason for 
this suppression of the current is that a fraction of the 
electron flow originating from the source is diverted into 
the tip due to the tip-wire coupling. While the condi- 
tion /t = ensures that the same electron current is 
re-injected into the wire, for symmetric tunneling the tip 
injects with equal probabilities right and left moving elec- 
trons. Hence half of the injected current flows back to 
the source electrode, causing the reduction of the two- 
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FIG. 3: (Color online) Zero temperature differential con- 
ductance as a function of the source-drain bias for a non- 
interacting wire with contact impurity strengths Ai = A2 = 
0.1 in presence of a tip with an applied voltage Vt adjusted to 
fulfill the condition /t = 0. Tunneling is totally asymmetric 
(x = 1) and the tunneling strength has the values 7 = 0.3 
(dotted line) and 7 = 0.7 (dashed line). The solid line repre- 
sents the case with 7 = 0. The result is independent of the 
tip position xq. The gate is grounded {Vo ~ 0), and the bias 
is applied symmetrically {Vs/o = ±^2). 

terminal conductance. As we shall see below, the situa- 
tion is different in the case of asymmetric tunneling. 

The second feature that can be observed in Fig. [2] is 
an alternating depth of the Fabry-Perot minima. This 
modulation originates from the interference between dif- 
ferent paths that arc possible for an electron ejected from 
the tip. For instance, the path of an electron ejected as 
right mover towards the drain can interfere with the path 
starting as left mover towards the source followed by an 
elastic backscattering at the source contact. The differ- 
ence in length between these paths corresponds to a new 
frequency in the oscillations, which causes the modula- 
tion of the peaks. In the case of Fig. [21 where the tip 
is located in the middle, this additional frequency equals 
twice the Fabry-Perot frequency, so that the tip affects 
every second minimum in the same way. As we shall see 
below, in general, the modulation pattern depends both 
on the asymmetry coefficient and on the position of the 
tip. The modulation effect arises to order 7^ A when we 
treat the impurity strength and tunneling amplitudes as 
perturbation parameters. 

The third effect of the tip consists in a reduction of 
the visibility of the the Fabry-Perot oscillations: in the 
presence of the tip the relative separation between max- 
ima and minima decreases. This reduction stems from 
the decoherencc introduced by the tip, since the proba- 
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FIG. 4: (Color online) Zero temperature differential conduc- 
tance as a function of the source-drain voltage for a non- 
interacting wire for several values of the tip position a;o = 
(solid line), xo — 0.17 (dotted line) and xq — 0.41 (dashed 
line). Tunneling has amplitude 7 = 1 and is symmet- 
ric (x = 0). The contact impurities have equal strengths 
Ai = A2 = 0.1. The gate is grounded (Vg = 0), and the bias 
is applied symmetrically (Vs/d ~ 



bility of constructive interference between paths with two 
backscattering processes at the contacts decreases when 
electrons can be incoherently absorbed and re-ejected by 
the tip. Notice that the reduction of visibility is of order 
7^A^, and it is therefore negligible with respect to the 
modulation effect in the Fabry-Perot regime. 

Let us now discuss the role of asymmetric tunneling in 
the voltage probe configuration. When x 7^ Oj the effect 
of conductance suppression is less pronounced then for 
symmetric tunneling. This can be seen already in the 
case of a clean wire (Ai = 0), where 

e2 2 + ^2 2 

iM^^iVs^Vu)^^^, (43) 
h 2 + 7 

and the value Vt of the tip voltage ensuring /t = is 
given by 

Vt = 1[Vs{1 + x) + Vb{1-x)] ■ (44) 

As one can see from the last factor in Eq. (j43|) . the 
suppression of the current Im is completely absent for 
fully asymmetric tunneling x = if- Importantly, this 
features persists also in the presence of realistic contacts 
(Ai 7^ 0), as shown in Figure [3l where the differential 
conductance dlyi/dV is plotted as a function of the 
source-drain voltage for several values of the tunneling 
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strength 7. Increasing the tunneUng strength simply de- 
creases the amphtude of the Fabry-Pcrot osciUations but 
does not change the average value of the conductance. 
Two more noteworthy features can be observed: in the 
fully asymmetric case also the modulation of the peaks is 
absent, and the non-linear conductance is independent of 
the tip position. The reason lies in the specific tunneling 
conditions. For example, a right moving electron ejected 
by the tip cannot be re-adsorbed after scattering as a 
left moving one, and this rules out interference effects 
between electrons traveling through the tip and electrons 
that have undergone an odd number of backscattering 
events at the contacts. Such processes would give rise 
to eflECcts related to the tip position, while interference 
phenomena with electrons that have undergone an even 
number of backscattering events, which continue to be 
present also for x — ^^rc independent of the tip 
position. Moreover, in the completely asymmetric case, 
electrons passing through the tip continue to move in 
the same direction, and this is the reason why, also for 
strong tunneling, the average value of the differential 
conductance is independent of 7. 

Finally, we analyze the dependence of the differential 
conductance on the tip position. For simplicity we limit 
this discussion to the case of symmetric tunneling illus- 
trated in Fig. m Apart from the conductance suppres- 
sion discussed above, one sees that the modulation effect 
exhibits a strong dependence on the tip position. In par- 
ticular, when the tip is close to a contact impurity, we 
observe Fabry-Perot-like oscillations over-imposed by an 
oscillation with large period due to coherent motion of 
carriers between the tip and the contact impurity remote 
from the tip. 



IV. THE INTERACTING CASE 

In this section we discuss the three terminal set-up 
in presence of electron-electron interaction. For arbi- 
trary values of the interaction strength, contact resis- 
tances, and tunneling amplitudes an analytical treatment 
is not possible, therefore we focus here on the Fabry- 
Perot regime. In this regime, characterized by highly 
transparent contacts to the electrodes, the role of in- 
teractions has so far only been analyzed for a two ter- 
minal set-up i^i^ Since the impurity strengths Ai are 
small, they can be treated perturbatively. The electron- 
electron interaction ([8]) will be accounted for exactly us- 
ing bosonization. The evaluation of the currents in the 
three terminals will be based on the out-of-equilibrium 
Keldysh formalismi^ We shall first discuss the effects of 
electron-electron interaction for the two-terminal set-up 
in the Fabry-Pcrot regime, i.e. in the absence of the tip, 
and then turn to the combined effect of tip and electronic 
correlations. 
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FIG. 5: (Color online) Differential conductance of the two- 
terminal set-up (in the absence of a tip) as a function of the 
source-drain voltage for several values of the interaction pa- 
rameter g = 1 (solid line), g = 0.75 (dashed line) and g = 0.25 
(dotted line). The contact impurities have equal strengths 
Ag 1 = Ag 2 ~ 0.1. The gate voltage is Vg = 0, and the bias 
is applied symmetrically (Vs/d = i'l^/S). 



A. Interaction effects on Fabry-Perot oscillations 
in a two-terminal set-up 

Let us first analyze the effects of electron-electron in- 
teraction for a contacted wire without tip. As in the non- 
interacting case, for 7± = the problem is reduced to a 
two-terminal set-up, where /t = and Ig = Id = Im- 
Furthermore, Im can again be written as a sum of the 
current Iq in a wire with adiabatic contacts and /imp, 
see Eq. (|26p . Importantly, while Iq is unaffected by the 
interaction in the wiref22i2^ the current Jimp, accounting 
for the contact resistances, is strongly modified by the 
interaction. One can still decompose /imp into 

eujt 

Amp = -rr- U'inc + jcoh) , (45) 
ZTT 

where jinc is the sum of two terms related to a single 
impurity each, and jcoh describes interference between 
scattering processes at the two impurities. Here, an 
important difference emerges with respect to the non- 
interacting case. The Fermi velocity in the wire is en- 
hanced by the interaction parameter g, leading to a 
higher ballistic frequency 



Moreover, the interaction is also affecting the strength of 
the contact impurities: the forward-scattering processes 
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FIG. 6: (Color online) The differential conductance dlm/dV [in units of (ea;2/27r)] of the two-terminal set-up (in the absence of 
a tip) is shown as a function of the dimensionless source-drain bieis u and the dimensionless gate voltage mg. The strengths of 
the contact impurities are equal, and are characterized by Af.i = 0.1 and Ae^i = 0.25. In panels (a), (b), and (c) the source-drain 
voltage is applied symmetrically, us/d = ±'u/2, and the interaction parameter is jr = 0.25,0.34, and 0.46, respectively. The 
dashed-dotted lines in panels (b) and (d) are a guide for the eye to identify the periodic pattern of the Fabry-Perot oscillations 
determined by the periods Ait and Auq- Their ratio yields the value of g, as shown in the table for the three cases. In panel (d) 
the source-drain voltage is applied asymmetrically (us = and ud = 0) to a wire with interaction strength g = 0.34. When 
compared with panel (b) for the same interaction strength, the asymmetric bias twists the pattern. 



are left unchanged whereas the backscattering ones are 
renormalized^ 

^^i>^h = ^<' . (47) 
[Ap.i = A 



where a-w — ayf/gL is a small dimensionless cutoff pa- 
rameter. The cutoff length aw, which is related to 
the lattice spacing or the electronic bandwidth of order 
fiuw/ow: is introduced in App. |B] 

In the Fabry-Pcrot regime we can again restrict our- 
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selves to terms up to third order in the contact impurity 
strengths A^. Then, the incoherent and the coherent con- 
tributions can be written as 



1=1,2 



and 



with 



and 

,■(2 



_ .(2) .(3) 
Jcoh — Jcoh + Jcoh ■ 



(48) 



(49) 



(50) 



Jcoh =2A^_iA^ 2^i2(u)cos{2[kw + .9(ww- wg)]} (51) 



Jcol =2A^.iA^.2(Af4+Af,2)Z?i2(w) 

X sin {2 [kw + 5("w - "g)]} , 



where we have introduced 



D,,{u) 



poo 

/ dr sin(uT) 
Jo 



(52) 



(53) 



sin[47rX**(e,;^,;r)] e^^^-o^ 



The dimensionless voltages u, uw and uq are now scaled 
by the factor e/hjuj}^ compared to the physical voltages 
— (Vs + T^d)/2 and Vg, respectively. In the ex- 
pression for Dij, the dimensionless integration time is 
defined as r = i^^t, and the functions 7?,*g(^; t) and 
X**(^;^';t) are the real and imaginary parts, respec- 
tively, of the auto-correlation function of the bosonic 
phase field $ introduced in App. |B] The quantity 
defined in Eq. (j53p is cut-off independent, since the cut- 
off dependence of the prefactor is compensated by the 
one of the correlation functions. Explicit results for the 
phase field auto-correlation function have been given in 
a previous paper iH Further, the £,i = Xi/ L {i = 1,2) are 
dimensionless contact impurity positions. Equations (|50f 
I52p are obtained from a perturbative development of the 
current in the impurity strengths A^ employing the meth- 
ods described in Apps. [X] and [B) The current j^^^ in 
Eq. ((5^ includes forward scattering processes that give 
rise to the factor Af,i + Af,2 and a twofold backscattering 
contribution leading to the factor Ag j^Ag 2- 

Another important effect of the interaction is that the 
incoherent term ji^c does not depend linearly on the 
bias as in the non-interacting case. Instead, it exhibits 
oscillations of period Au ~ tt, due to the interplay 
between backward scattering at one contact impurity 
and Andreev-type reflection at the other contact On 
the other hand, the coherent term jcoh, responsible for 
Fabry-Perot oscillations, shows a power-law suppression 
with increasing voltage . ^^'^^ Thus, in the presence 
of interaction two types of oscillations are present, 
namely the Fabry-Pcrot ones (already existing for a non- 
interacting wire and modified by the interaction), and 



the Andreev-type ones (purely due to the interaction). 
These two types of oscillations are characterized by the 
same period in the source-drain bias, and they are of 
the same order in the impurity strength, if we assume 
that the two contact transparencies are comparable 
(Ai ~ A2). It is therefore difficult to distinguish the 
two phenomena from an inspection of the two-terminal 
differential conductance, which is shown in Fig. O as a 
function of the source-drain bias for various values of the 
interaction parameter g. Besides the power-law suppres- 
sion of the amplitude at high applied bias, we see that 
for strong interaction {g < 1/2) the sinusoidal behavior 
of the oscillations is deformed into a saw-tooth-like 
shape. Furthermore, although the total current (pS)) in 
the presence of contact resistances is always smaller than 
the current Iq of an ideally contacted wire (/imp < 0), 
the differential conductance may exceed e^//i. This is a 
well known effect of non-linear transport in Luttinger 
liquids^, reflecting the fact that the conductance cannot 
be expressed in terms of single-particle transmission 
coefficients. In Sec.|V]we shall comment on how the two 
types of oscillations may be experimentally distinguished 
in a three-terminal set-up. 



Further interesting insights emerge from the analysis 
of the conductance dlu/dV as a function of both the 
source-drain bias V = Vs — Vu = {huj1/e)u and the 
gate bias Vg = {tiuj*j^/ e)uQ. Corresponding conductance 
plots are shown in Fig. [B] Panels (a), (b) and (c) refer 
to three different values of the interaction strength g, in 
the case of a symmetrically applied source-drain bias, 
ug/D = ±u/2. The oscillations of the conductance as a 
function of V and Vg are characterized by two periods 
AV^ and AVg- The period Ay coincides with the 
period of the function Di2{u) [Eq. (|53|l ] appearing in the 
coherent terms ((5T|) and ((52|) . since the functions Dii(u) 
and D22(u) related to the incoherent contribution ([50]) 
exhibit the period Ay/2. We thus recover the result 
of Ref. |2^. On the other hand, the period AVq in the 
gate voltage is determined by the sinusoidal factors 
of Eqs. dHI) and The values of AV and AVb 

depend on the interaction strength g and are inversely 
proportional to g and g^, respectively. Interestingly, 
the ratio of these periods yields the Luttinger liquid 
interaction strength, ATZ/AVg = Am /Aug = 2g, as can 
be checked from the table associated with Fig. [B] 
Panel (d) describes the case of an asymmetrically applied 
bias (ug ~ u and ud = 0), for the same interaction 
strength as panel (b). In this case uw = u/2 [see 
Eg. ([5^ ]. so that an additional dependence on V arises 
from the sinusoidal factors of Eqs. ((5T|) and ([52]) . and 
the period in V at fixed Vq changes. For this reason the 
two-dimensional pattern of the nonlinear conductance 
is twisted with respect to panel (b). However, the 
quantities AV and AVg related to a symmetrically 
applied bias can still be obtained, 6.17., by projecting 
the conductance maxima on the y-axis and measuring 
the distance between these projections as indicated by 
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the arrows in panel (d). The value of g can therefore be 
extracted also in this case as Am /Auq = 2g. We remark 
that a qualitatively similar twist of the conductance pat- 
tern has recently been observed in carbon nanotubesJ^ 

Conductance plots as a function of the transport 
and gate voltages have previously been discussed in 
the context of carbon nanotubes in Rcfs. 2il and [2^. 
We point out that the way we introduce the bias and 
gate voltages in our model [see Eq. differs from 
the one adopted in the above papers. Our approach 
accounts for several basic physical facts. In a non-chiral 
quantum wire only the electrochemical potentials of 
the leads can be controlled experimentally, whereas the 
electrochemical potentials of right and left movers inside 
the wire are a result of the biasing of the wire and its 
screening properties. As a consequence, the source and 
drain biases, Vs and Vd, are applied here only in the 
related leads. This is in accord with a basic hypothesis 
underlying the definition of an electrode, namely that 
inelastic processes in the lead equilibrate absorbed 
electrons, yielding a voltage drop at the contacts even in 
the absence of contact impurities. On the other hand, 
the charge density of metallic electrodes is typically 
insensitive to a gate, due to their electroneutrality. For 
this reason, in our model the gate voltage Vq is applied 
only to the interacting wire and not to the leads. 
The precise form of the coupling to the biasing voltages 
adopted in the model has implications on the behavior 
of the current as a function of bias and gate voltages. 
We find that the dependence on Vq and (Vs + ^d)/2 
involves a factor g^, as shown, for instance, in Eqs. (j51[) 
and ([5^ . [In the dimensionless formulation one factor 
of g is contained in the definition of the dimensionless 
quantities uq and uvv-] The difference (Vs + Vd)/2 — Vq 
is proportional to the bare electron charge injected 
into the wire, whereas the g^ factor originates from 
the partial screening occurring in a Luttinger liquid)^ 
and physically describes the fraction of the bare charge 
that remains unscreened. In particular, in the limit 
(yi — > of an electroneutral wire we obtain that the 
current depends only on the difference Vs ~ and is 
independent of the gate, as it should be. 

On a more formal level, these physical properties are 
encoded in the zero modes '^o,±{x) [see Eq. (jB13|) ]. In- 
deed, the transformation $± — > (f>± 4- <f>o.± of the chiral 
boson fields gauges away the bias term ([5]). We note that, 
differently from the homogenous Luttinger liquid case, 
in the presence of leads the zero modes cannot be just 
linear functions of the position uniformly along the en- 
tire system. The inhomogeneity of the system leads to a 
non-trivial space dependence of the zero modes '^o,±{x), 
which can be obtained from the boson Green function of 



the inhomogeneous LL model, as shown in Eq. (jB14p . 
B. Interaction effects on electron tunneling from 
the tip: The case of adiabatic contacts 

We shall now consider the full three-terminal set-up, 
and discuss the effects of the wire electron-electron in- 
teraction on tunneling from the tip, both for the case of 
electron injection and in the voltage probe configuration. 
We start by presenting results for a wire with adiabatic 
contacts (Ai = 0). For a non- interacting wire, the calcu- 
lation described in Sec. IIIII vields 

^^^7r(2TWr^"^ '"^^J (54) 
and 

= ^(^^^{l^s [4 + 27^(1 -X + 7V)] (55) 

-Vt,[A + 27^(1 + X + 7'x')] + 4x7'^t} , 

where 7 is the total tunneling strength defined in Eq. ([34]) 
and X is the tunneling asymmetry parameter introduced 
in Eq. (j35p . Thus, in the absence of interaction, the 
currents depend linearly on the three applied voltages 
and are independent of the position of the tip. 

When electron-electron interaction is taken into ac- 
count, an exact solution of the tunneling problem is 
not possible for arbitrary values of the tunneling ampli- 
tudes 7±. We shall assume that 7± ^ 1, consistent with 
the tunnel Hamiltonian approach, and provide results to 
leading order in perturbation theory. The currents in the 
source and drain leads are again written as in Eqs. (jl6p 
and (fT7|) , where /m and It are evaluated now to order 7^ 
yielding 

= I0+ hi.-f^ (56) 

and 

It = /t,7^ (57) 

where 

Im{T),^^ = -^(7*)^ JM(T),72 ■ (58) 

Here 

8+8-1-2 

7 =7 aw " (59) 

is the tunneling amplitude renormalized by the electron- 
electron interaction. The dimensionless currents jM{T}.-y^ 
read 
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where 



^ /drgM(T)(r)sin{47r[l*-*+(6;Co;r)+X^nO;0;r)]}e^^L^- «o;Co;.)+7^-(0;0;.)J ^^^^ 

5 Jo 



Qm(t) — sin(MT/2) COS [(uw — ut)t] + X cos(ur/2) sin [(uw — ut)t] 

Qt{t) = 2cos(mt/2) sin [(ww — wt)t] + 2xsin(ur/2) cos [(ww — 'Ut)t] 

I 



(61) 



Here is a small dimcnsionlcss cutoff parameter for 
the tip defined in App. [B) The functions 7?.reg*^ {CiC'i t) 
and 2r*+*+ (^; r) arc the real and imaginary parts 
of the auto-correlation function of the chiral wire field 
defined in Eqs. (|C7p and (jC8|) . respectively, while 
^reg(C;^V) and X'^'^(^;^';t) are the real and imagi- 
nary parts of the correlator of the tip field Lp given in 
Eqs. jCTTI) and (|CT2|) . The integral ^ is a cut-off in- 
dependent quantity. 

We consider two parameter domains of the three- 
terminal set-up corresponding to the cases where the tip 
operates as an electron injector and as voltage probe, 
respectively. In the electron injection case, source and 
drain are at the same electrochemical potential while a 
bias is applied to the tip. For this configuration the cur- 
rent noise was evaluated in Refs. [l3 and|l^. Here we shall 
explicitly evaluate the non-linear tunneling conductances 



equality 
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DT 
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(63) 
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Conventional Luttinger liquid theory, where the pres- 
ence of the source and drain electrodes is neglected, pre- 
dicts that an electron charge injected by tunneling e.g. as 
a right-mover into an interacting wire breaks up into sep- 
arate charge pulses moving in opposite directions, namely 
a fraction (1 + g)/2 moving to the right and a fraction 
(1 - g)/2 going to the left J^i^i25i36,37^ r^j^.^ ^^^^^ 

inates from the coupling between the densities of right 
and left moving electrons, accounted for by the homoge- 
neous LL Hamiltonian. As a consequence, one expects 
that when the tip injects electrons asymmetrically, e.g. 
only toward the drain electrode on the right (x = 1), the 
electron-electron interaction would cause a part of the 
current to flow also to the source electrode on the left. 

However, when the source and drain electrodes are ex- 
plicitly taken into account, our results show that the 
above expectation is in fact wrong. Remarkably, using 
Eq. ([50]) . one can indeed prove that for Vs = Vd the 



Aix) 



l/s 



I^d| + |/sI 



X 



(64) 



holds, indicating that for a clean wire the current asym- 
metry is independent of the wire interaction strength g. 
In particular, for fully asymmetric tunneling (x = 1), 
the whole current is injected into the drain electrode, 
just as in the non-interacting case. This unidirectional 
charge flow even in the presence of interaction arises from 
the phenomenon of Andreev-type reflections Even 
though charge fractionalization occurs in the bulk of 
the wire, the plasmonic excitations reaching an interface 
with the leads experience the mismatch of the interaction 
strengths in the wire and in the electrode and are thus 
partly reflected as an oppositely charged excitation. The 
sum of all reflected pulses at both interfaces restores the 
property that the whole current flows into the drain, like 
in the non-interacting wire. This behavior is in fact very 
similar to an effect occurring in a two-terminal set-up, 
where the conductance of a wire adiabatically connected 
to electrodes is G2t = independent of the interac- 

tion strength. Thus, for perfectly transmitting contacts, 
it is impossible to extract the interaction constant neither 
from the conductance of a two-terminal set-up nor from 
the current asymmetry in three-terminal measurements. 

Nevertheless, in a three-terminal set-up signatures of 
interaction do appear in the behavior of the differential 
conductances Gst and Gdt as a function of the tip- 
source and tip-drain bias. Figures [Zla) and [Tfb) show 
Gst and Gdt for the case of a tip located in the middle 
of a wire with interaction strength g = 0.25. The var- 
ious curves correspond to different values of the asym- 
metry parameter x, which unbalances the amount of in- 
jected right vs. left moving electrons. The fully sym- 
metric case (x = 0) was discussed in Ref. H^. While 
for a non-interacting wire Gst and Gdt are constant [as 
can easily be seen from Eqs. ((54)) and (|55|) ]. in the pres- 
ence of interaction an oscillatory behavior arises. These 
oscillations are entirely due to the electron-electron in- 
teraction in the wire, which causes Andreev-type reflec- 
tions even at adiabatic contacts. With increasing x the 
conductance Gst decreases until it vanishes for x = Ij 
whereas the conductance Gdt increases up to the max- 
imum value for the completely asymmetric case. The 
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FIG. 7: (Color online) Electron injection into an interacting wire. Panel (a) [(b)] shows the tunneling diflerential conductance 
Gts [Gtd] between the tip and source [drain] electrode as a function of the dimensionless tip-source [tip-drain] bias for a wire 
with interaction strength g — 0.25. The various curves refer to different values of the tunneling asymmetry (x = solid line, 
X = 0.6 dotted line, and x = 1 dashed line). The tunneling strength is 7* = 0.01, and the tip is located in the middle of the 
wire. Panels (c) and (d) are the same as panels (a) and (b) but the tip is located near a contact at xq — 0.451/. 



relation Gst = Gdt(1 — x)/(1 + x) between these two 
conductances is independent of g. 

Figures [7]^c) and[71Jd) describe the case of an off-centered 
tip located at xq = 0.45L. Apparently, the period of the 
oscillations is the same as in panels (a) and (b) where 
the tip is in the middle. This is due to the fact that this 
period is related to the traversal time of plasmonic ex- 
citations originating from the tip and interfering at the 
same point after an even number of Andreev-type reflec- 
tions at the contacts. This traversal time depends neither 
on xo nor on the asymmetry coefficient. 

Let us now discuss the configuration where the tip acts 
as a voltage probe, i.e. when Vs 7^ Vu and Vr is set to a 
value such that /t = is fulfilled. In this configuration 
the quantity on the left hand side of Eq. (j36p is vanishing, 
due to Eqs. (fT6|) and (fT7|) . By applying a source-drain 
bias, one can analyze how the source-drain conductance 
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SD 



(65) 



is affected by the interaction strength. It is worth em- 
phasizing that in a two-terminal set-up, i.e. in the ab- 
sence of the tip (7± = 0), one obtains for a clean wire 
GsD ~ G2t ~ e^jh , independent of the interaction 
strength. As already mentioned previously, this is due 
to the fact that, although the electron charge injected by 
the source splits up in fractions through the interaction- 
induced Andreev-type reflections at the contacts, in a 
clean wire the series of these fractions always sums up 
to e, disguising the interaction effects in the dc average 
current Our results show that a quite different behav- 



ior emerges for a three-terminal set-up, even in the con- 
figuration where the tip does not inject any net current 
into the wire. Figure[5I^a) shows Gsd as a function of the 
source-drain bias, for different values of the interaction 
strength, ranging from a non-interacting to a strongly 
interacting wire. The left panel refers to the case of 
symmetric tunneling x — whereas the right one an- 
alyzes the role of a tunneling asymmetry. As one can 
see, the effects of interaction in the wire become observ- 
able through the voltage probe, since oscillation of Gsd 
originating from Andreev-type reflections emerge. Notice 
that at constant bare tunneling amplitude 7 the zero bias 
conductance is higher in the presence of interaction than 
for a non-interacting wire, since the renormalization ((59|) 
of the tunneling amplitude suppresses 7. With increasing 
tunneling asymmetry [sec Fig. El^b)], the differences be- 
tween interacting and non-interacting wires become less 
pronounced, and indeed the oscillations are washed out 
for fully asymmetric tunneling x = il- 

The dimensionless tip voltage mt = (e/fiw^)!^ 
ensuring /t=0 shows an interesting dependence on the 
source-drain bias. In the limiting cases of symmetric and 
completely asymmetric tunneling this dependence coin- 
cides for interacting and non-interacting wires (namely 
ut ~ uw for X = and ut = mw i u/2 for x = il)- 
For intermediate values of the asymmetry parameter x 
the tip voltage ut shows an oscillatory behavior with 
period Aut = 27r as a function of the source-drain bias. 
We also see that the period of Gsd in Fig. [5] is twice as 
large as the period of Gts and Gtd in Fig. [7] where the 
tip is in the electron injection configuration. This is due 
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FIG. 8: (Color online) The effect of the tip in the voltage 
probe configuration on the zero temperature conductance is 
shown as a function of the source-drain bias. Panel (a): the 
case of symmetric tunneling (x = 0). Panel (b): the case of an 
asymmetry in tunneling (x = 1/2). Different curves refer to 
different values of the interacting strength: non-interacting 
{g — 1, solid curve), weakly interacting {g = 0.7, dashed 
curve), moderately interacting {g = 0.4, dotted curve), and 
strongly interacting {g = 0.25, dashed-dotted curve). The 
bare tunneling strength is 7 = 0.5 and the dimensionless cut- 
off parameter is aw ~ 10^''. 



to the fact that in Fig. [5] the source-drain bias is applied 
symmetrically (uw — 0) while in Fig. [7] source and drain 
are both grounded and the bias is only applied to the 
tip. 

Finally, we emphasize again the difference between the 
electron injection and the voltage probe configurations 
of the tip: While in the former case an asymmetry 
in tunneling does not spoil the observation of effects 
of electron-electron interaction (sec Fig. [7|), in the 
latter case interaction-induced oscillations can be best 
observed for symmetric tunneling and they are in fact 
vanishing for fully asymmetric tunneling. 



tion, contact impurity scattering and electron tunneling 
from the tip. In particular, we discuss how a finite con- 
tact resistance modifies the Andreev-type oscillations of 
the tunneling conductances, previously discussed for the 
case of adiabatic contacts (A^ = 0). We present results 
obtained by perturbation theory for weak contact impu- 
rities Ai, and tunneling amplitudes 7±. Technical details 
can be found in the Appendices. The currents may be 
written as 



'M,7 



2 + /] 



M,72A 



(66) 



and 



(67) 



Here Iq = (e^//i)(Vs — Vd) is the current of an ideally 
contacted wire in the absence of the tip, whereas /imp 
is the leading order term accounting for non-ideal con- 
tacts [see Eq. (|45|) ]. In general, this latter term involves 
both Fabry-Pcrot and Andreev-type oscillations. Both, 
Iq and /imp, vanish when the electrochemical potentials 
for source and drain electrodes are equal (Vs = Vd); al- 
ternatively, they can be easily determined by measuring 
the current-voltage characteristics in the absence of the 
tip. Henceforth, we shall focus on contributions to the 
currents arising from the presence of the tip. The lead- 



ing order terms Im,^^ and /t,72, given by Eq. ([58|) . de- 
scribe tunneling into an ideally contacted wire and con- 
tain only Andreev-type oscillations. The next-to-leading 
order terms (7^ A, 7^A^, ...) also exhibit oscillations orig- 
inating from interference between backscattering at the 
contacts and tunneling to/from the tip. Such oscillations, 
though modified by the interaction, are already present 
in a non-interacting wire, unlike the Andreev-type oscil- 
lations of the leading order terms (7^), that are instead 
entirely due to the interaction. We thus analyze how in- 
teraction affects the terms /t(m),72Aj which represent the 
most relevant correction to the Andreev-type oscillations 
discussed above. These terms describe to leading order 
the interplay between electron injection at the tip and 
backscattering at the S and D contacts. 



Explicitly one finds 



C. Interaction effects on electron tunneling from 

the tip: The case of a wire with non-ideal contacts /m(T) -y^X — — ^(7*)^ A* •\/l — JM(T) j^X (68) 

In this section we analyze the three-terminal transport 

properties in the presence of electron-electron interac- where A* = Ag + Ag 21 a-i^d 
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FIG. 9: (Color online) Zero temperature non-linear tip-drain conductance Gdt as a function of the tip bias for tunneling 
amplitude 7* = 10~^ and symmetric tunneling x ~ 0. The upper panels (a), (b) and (c), are related to a tip in the middle of 
the wire (xq — 0), whereas the lower ones (d), (e) and (f) to a tip located at xq ~ 0.451/. Panels pairs (a) and (d), (b) and 
(e), and (c) and (f) describe the case of a wire with weak {g — 0.75), moderate {g — 0.5), and strong {g — 0.35) interaction 
strength, respectively. In each panel the different curves refer to different contact impurity strengths. The solid curves describe 
the case of ideal contacts where the oscillations are purely Andreev-type. The dashed [dotted] lines refer to finite contact 
impurity strength Ab,i = Xb 2 = 0.1 [= 0.2]. For a weakly interacting wire the conductance oscillations are mostly due to 
the conventional interference between backscattering at the contacts and tip tunneling, and Andreev-type oscillations become 
visible only for extremely low contact resistance. In contrast, for stronger interaction strength a finite contact resistance is 
sufficient for the oscillations to be attributed to Andreev-type processes. The inset of panel (e) shows the definition of the 
average amplitude referred to in the text. 
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(69) 



r 



where 

=V2+Vi- 2viV2 (70) 

pr^ = 2(,^2-m). (71) 

The functions F^^^f^ (n; T2) and F^^^Jj', ^{n; T2) are 



defined in App. [B]in Eqs. (|B8)) and (|B24p . respectively, 
and the functions 7?.w(Co; Ci; n; 7-2), Iw(Co; n; 7-2), 
TZt{ti; T2) and Xr(Ti; T2), accounting for the real and the 
imaginary parts of several correlation functions in the 
wire and in the tip, are defined in App. \C\m Eqs. (|C1[) . 
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((C2|) . (|CTT|) and (|CT2|) . respectively. 

For simplicity, we limit the analysis of Eqs. (j68p and 
(j69p to the electron injection configuration where source 
and drain are grounded. We start with the case of sym- 
metric tunneling (% = 0). 

As already observed in Sec. IIIII for a noninteracting 
wire, the term ([55]) leads to additional oscillations in 
the differential conductance of the three-terminal set-up. 
These conventional oscillations are characterized by two 
periods related to the distances between the tip and 
the contact impurities, so that the pattern depends on 
the tip position. Electron-electron interaction modifies 
this pattern reducing the amplitude of the conventional 
oscillations and giving rise to additional Andreev-type 
oscillations. The case of a tip in the middle of the wire 
is shown in the upper panels (a), (b) and (c) of Fig. [9l 
where the differential conductance Gdt is plotted as a 
function of the tip bias ux for three different values of 
interaction strength, ranging from weak (17 = 0.75), over 
moderate (g = 0.5) to strong interaction {g = 0.35), 
as displayed in the three panels. In each panel the 
solid curve refers to the case of ideal contacts where the 
oscillations are purely of Andrcev-typc. The dotted and 
dashed curves describe the effect of finite contact resis- 
tances arising from the contribution of the term ([S5)) . 
As one can see from panel (a), for weak electron-electron 
interaction the conventional oscillations dominate and 
mask the Andreev-type oscillations. In this case, only 
extremely good contacting might allow to identify 
Andreev-type processes. However, for moderate inter- 
action strength [panel (b)], the two types of oscillations 
have comparable amplitudes, and for strong interaction 
[panel (c)] the conventional oscillations are strongly 
suppressed while the term ((55]) only causes a small shift 
of the conductance value. The oscillations of Gdt are 
essentially Andreev-type. 

A similar effects occurs when the tip is closer to one 
of the contacts displayed in the lower panels (d), (e) 
and (f) of Fig. [S) The main difference is that in this 
case the pattern of the Andreev-type oscillations is more 
sinusoidal, even for weak interactions. 

Our result indicates that, for a wire with a given in- 
teraction strength, there is crossover value of the 
(renormalized) contact resistance, below which the os- 
cillations of the non-linear conductance can essentially 
be attributed to Andreev-type processes. We have quan- 
tified Aq for the case of a tip close to the contacts, where 
the regularity of oscillations allows for a straightforward 
determination of their amplitude, defined as the average 
distance between maxima and minima, as schematically 
displayed in the inset of Fig.[SJe). The crossover impurity 
strength is then simply determined by the value of A* 
for which the amplitude Ary2x of the conventional oscilla- 
tion term lT,'y^x [see Eq. (|68p ] equals the amplitude A^2 
of the Andreev-type oscillation term It^j2 [see Eq. ([55]) ]. 
The result is given in Table U for different values of inter- 
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0.5 


0.6 


0.7 


0.8 


Ac 


0.2 


0.02 


4 • lO"^' 


4 ■ 10-* 


10-* 



TABLE I: Crossover value of the (renormalized) contact im- 
purity strength Ac, below which oscillations can be attributed 
to Andreev-type processes, for various values of the interac- 
tion strength g. The tip is located at xq — 0.45L, like in the 
lower panels of Fig[9] 

action strength. For contact impurity strength A* < Xq 
the oscillations of the non-linear conductance are essen- 
tially of Andreev-type. 

Let us finally briefly consider the case of asymmetric 
tunneling % 0. An important result is that, in view 
of Eq. ([55]) . the contribution to the current of order 7^ A 
vanishes in the case of totally asymmetric tunneling (x = 
±1). This property is thus robust to electron-electron in- 
teraction within the Luttinger liquid picture. In fact, one 
can show that in this case only perturbative contributions 
of order 72"(Ai A2)"+™ (n = 1, 2, 3, . . . ; m = 0, 1, 2, . . .) 
are nonvanishing. 

V. DISCUSSION AND CONCLUSIONS 

In the present work, we have investigated transport 
properties of a quantum wire contacted to source and 
drain reservoirs in the presence of a third electrode 
(tip) injecting electrons into the wire. We have tailored 
our model to account for various aspects of a typical 
experimental situation by including finite contact 
resistances and the presence of a gate in addition to 
electron-electron interaction, and by analyzing the 
effect of the position of the tip as well as the role of a 
tunneling asymmetry. Specifically, we have considered 
both the situation where the tip behaves as an electron 
injector and the voltage probe configuration. We have 
found that the three-terminal set-up exhibits extremely 
rich behaviors, determined not only by each of the 
above aspects, but also by their interplay. In order 
to facilitate the discussion, we propose to the reader 
different perspectives from which our results can be 
considered. 

The effects of electron- electron interaction on Fabry- 
Perot oscillations. The origin of Fabry-Perot oscillations 
boils down to quantum interference between electron 
backscattering at two (or more) impurities. As a 
consequence, this phenomenon is present also in a 
non-interacting quantum wire (see Sec. IIIip . where the 
oscillations appear both as a function of the source-drain 
bias and as a function of the gate voltage. The interfer- 
ence pattern is modified by electron-electron interaction, 
which introduces a power-law suppression of the ampli- 
tude and, especially for g < 1/2, deforms the sinusoidal 
shape towards a saw-tooth-like shape (see Fig. 
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Interaction also leads to a (partial) screening of the 
charge in the wire^i, causing a change of the oscillation 
period as a function of the gate bias with respect to the 
period as a function of the source-drain bias. This effect 
suggests an operative procedure to extract the Luttinger 
liquid parameter g from measurements of the non-linear 
conductance in the Fabry-Perot regime (see Fig. [S]). The 
effects of an asymmetrically applied source-drain bias 
have also been discussed. We emphasize that, differently 
from previous approaches adopted in the literature, our 
way to introduce the biasing voltages correctly recovers 
both gauge invariance^ and the property that, in the 
limit of strong interaction t/ — > 0, the current-voltage 
characteristics only depends on the difference between 
source and drain bias Vs — Vd. 

Conventional vs. Andreev-type oscillations. Besides mod- 
ifying Fabry-Perot oscillations, electron-electron interac- 
tion also yields another major effect, which is absent 
in a non-interacting wire: At the wire-electrode inter- 
faces, plasmon excitations arc partially reflected due to 
the mismatch of interaction strengths in the interacting 
wire and the non-interacting electrodes. This effect, en- 
tirely due to interaction, occurs also for ideally contacted 
adiabatic interfaces and gives rise to a different type of os- 
cillations, which are termed Andreev-type oscillations^ 
since the incoming charge and the fractional charge re- 
flected at the contact have opposite signs, just as at an 
interface between a normal metal and a superconductor. 
In real experiments with interacting quantum wires in 
the Fabry-Perot regime, the current-voltage characteris- 
tics will in general exhibit both conventional Fabry-Perot 
oscillations, i.e. oscillations that are already present in 
a non-interacting wire and that are simply modified by 
interaction, and Andreev-type oscillations, purely origi- 
nating from interaction. The interesting question arises 
whether one can distinguish between these two oscilla- 
tory phenomena in an operative way and, in particular, 
whether it is possible to determine regimes and condi- 
tions, under which the latter can be observed. 

Since the amplitude of Fabry-Perot oscillations is 
roughly proportional to the reflection coefficients of the 
contacts whereas Andreev-type processes occur even with 
ideal interfaces, one might at first think that with improv- 
ing transparency of the contacts the non-linear conduc- 
tance of a two-terminal set-up would exhibit a predomi- 
nance of Andreev-type oscillations over the conventional 
Fabry-Perot ones. This is, however, not the case, since 
for an ideally contacted wire the sum of all Andreev- 
type refiection processes at the two interfaces exactly re- 
covers the injected pulse, when the sign of all reflected 
charge pulses is taken into account. The transmission 
of an interacting wire adiabatically connected to non- 
interacting leads turns out to equal 1, as was pointed 
out in Refs. [13 and [H. Although Andreev-type oscil- 
lations of the conductance do appear in the presence of 
even a single impurity^, their amplitude is proportional 
to the impurity reflection coefficient. This implies that 



two-terminal set-ups are not suitable to distinguish be- 
tween Andreev-type and Fabry-Perot oscillations, since 
both oscillations have the same dependence on the im- 
purity strengths Ag ^ . Furthermore they also exhibit the 
same period as a function of the source-drain bias. 

In contrast, our analysis suggests that a three-terminal 
set-up may allow one to distinguish Andreev-type oscilla- 
tions from conventional oscillations. As far as Andreev- 
type oscillations are concerned, three-terminal set-ups 
indeed offer one important advantage with respect to 
two-terminal ones: in the presence of a third electrode, 
Andreev-type oscillations appear even for the ideal case 
of a wire adiabatically connected to the source and drain 
electrodes (Ai=0). In the presence of interaction the tip- 
source and tip-drain non-linear conductances Gst and 
Gdt oscillate as a function of the tip voltage ut already 
to leading order 7^ in the tunneling amplitude, indepen- 
dent of contact impurity strengths A^. This effects holds 
when the tip acts as an electron injector (sec Fig. [T]) 
as well as when it acts as a voltage probe (see Fig. [8|) , 
and the oscillations vanish for a non-interacting wire [see 
Eqs. ([51]) and (|55|) ]. Thus, quite differently from a two- 
terminal set-up, in three-terminal set-ups Andreev-type 
oscillations become better visible when the contact trans- 
parency is improved. 

In view of the fact that in realistic experiments the 
contact resistance is always finite, we have quantitatively 
evaluated the infiuence of the contact resistance on the 
conductance oscillations [see Ea. ([S9)) ] showing that 
additional Fabry-Perot-type oscillations superimpose 
with the Andreev-type ones (see FigE]). We have thus 
put forward criteria for observing the interaction induced 
Andreev-type oscillations. At least two experimental 
situations are promising: For the conventional case of 
symmetric tunneling form the tip, we have determined 
typical values of the contact resistance below which the 
oscillations in the current-voltage characteristics can 
essentially be attributed to Andreev-type phenomena. 
The result, shown in Table [H indicates that the stronger 
the interaction of the wire the larger are the contact 
resistances that are tolerable to still observe Andreev- 
type oscillations. Furthermore, in case that the set-up 
allows for fully asymmetric tunneling, the leading order 
correction (|69p competing with the Andreev-type term 
is vanishing, even in the presence of interaction. 
In summary, in systems like carbon nanotubes where the 
interaction strength is typically strong, g ~ 0.2 — 0.3, 
while electron injection from an STM tip is typically 
symmetric, Andreev-type oscillations may be observed 
by achieving a high quality of the contacts to the leads. 
In contrast, in semiconductor quantum wires, where the 
interaction strength is usually moderate g ~ 0.6 — 0.7, 
asymmetric tunneling induced by a magnetic field is 
more suitable to observe Andreev-type oscillations. 

The effects of asymmetric tunneling. The above- 
mentioned case of asymmetric tunneling deserves some 
further remarks. Recent experiments by Yacoby and co- 
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workersiSiSi have shown that fully asymmetric tunneling 
into semiconductor-based quantum wires can be realized 
by appropriate tuning of a magnetic field. Inspired by 
these experiments, we have considered the possibility of 
an asymmetry in electron tunneling from the tip. Be- 
fore discussing our results we would like to point out the 
relation between our model and Yacoby's experimental 
set-up. While Yacoby et al. study electron tunneling 
between two parallel wires where momentum conserva- 
tion is required, our model considers injection from a 
point-like tip. Although these two situations may at first 
seem incompatible, a regime can be determined where 
they are equivalent. In the experiments of Refs. IT3l2ll 
electrons are injected from an upper shorter wire with 
length Lu into a lower longer wire with length Li. Since 
the tunneling region reasonably coincides with the length 
of the short wire, momentum conservation only holds up 
to an uncertainty 5k ~ l/L^- Although this uncertainty 
is small enough to select a specific electron momentum 
state in the upper wire, 5k may be much bigger than the 
mean level spacing of the lower wire, if the latter is much 
longer than the former (L/ 3> iu)- In this regime, while 
the electron wave function behaves like a plane wave for 
the short wire, for the long wire it can effectively be con- 
sidered as a localized wave packet, and our model applies. 

Under these conditions several interesting effects 
emerge. In the first instance, by using the tip as an elec- 
tron injector, the tunneling asymmetry can be exploited 
to gain the transmission coefficient of each contact by 
measuring the current asymmetry (|36p in the two cases 
of tunneling purely to the right (x = +1) and to the left 
(x = — 1), as has been shown in Eq. (HO]). Secondly, when 
the tip is used in the configuration of a voltage probe, 
fully asymmetric tunnelling allows to eliminate the sup- 
pression of the source-drain conductance Gsd , which oc- 
curs for symmetric tunneling. Similarly, Gsd becomes 
independent of the tip position. 

When electron-electron interaction is taken into ac- 
count, the scenario is even richer. Luttinger liquid the- 
ory predicts that electron-electron interaction induces 
a current asymmetry which depends on the interaction 
strength g. The appealing question arises whether this ef- 
fect is observable in experiments, where currents are mea- 
sured not directly in the interacting wire but in metallic 
electrodes connected to it. The investigation carried out 
in Ref. [13, based on the assumption that the interfaces 
between the interacting wire and the electrodes can be 
treated phenomenologically with a transmission coeffi- 
cient a la Landaucr-Biittiker, has led these authors to the 
claim that the interaction strength can be observed via 
the current asymmetry. Here we have scrutinized this 
prediction by taking the presence of source and drain 
electrodes into account fully consistently within the in- 
homogeneous Luttinger liquid model. Considering as a 
test bench the case of a wire adiabatically contacted to 
source and drain electrodes, we have proven that, al- 
though charge fractionalization does occur in the bulk 
of the wire, the sum of Andrcev-type reflection processes 



at the contacts leads to a current asymmetry A that is 
independent of the electron-electron interaction strength, 
just as it is the case with the two terminal conductance 
G-zt- Thus, already for this ideal case, no proof of charge 
fractionalization can be gained from the analysis of A, or 
from the ration e'^A/{hG2t)- We have also shown that, 
nevertheless, interaction effects do appear in the behavior 
of the nonlinear conductance, where interaction induced 
oscillations arise as a function of the tip-source and tip- 
drain bias. It is worth emphasizing that this feature is 
due to the three-terminal set-up, since the two-terminal 
conductance of a Luttinger liquid ideally contacted to 
leads is independent of the source-drain bias. 
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APPENDIX A: KELDYSH FORMALISM AND 
PERTURBATIVE EVALUATION OF THE 
CURRENT 

In order to compute the current in the three terminal 
set-up, we adopt the Keldysh formalism^ suitable to 
account for out-of-equilibriuni properties. According to 
Eq. (fTi|) . the current at position x (located in the source 
or in drain leads) and time t can be written as 

r;=± 

where 

j'^"\x,t) = J2r:-^l^'J\x,t)-^i''\x,t) : . (A2) 

Here f] = + {r] = —) corresponds to the upper (lower) 
branch of the Keldysh contour depicted in Fig. [TOl The 
current I{x,t) and various other quantities introduced 
below also depend on the injection point xq and the impu- 
rity positions xi and X2. These variables will frequently 
be suppressed to simplify notation. 

+ 



i > 



FIG. 10: Keldysh contour along the time axis. 
In the Keldysh interaction picture with 

Ti-O = Tlkin.yV + "Hu + "^kin.T (A3) 
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and 

one obtains 
I{x,t) 



'Hi — H\ + Htun + + H.fj.rj. , 



(A4) 



(A5) 



J2 (tk [j(")(x,t)c-*^./^±'''/rL^*'«<''''(ol\ 

r,=± \ L -I / 



is the Kcldysh time-ordering operator. Expanding the 
exponent in Eq. (|A5|1 perturbatively in terms of 7 and 
A, one obtains the eurrent to the desired order. Below 
we sketch the calculation of /^2_\, i.e., the contribution 
of order 7^ A to I{x, t). With the abbreviations 



Oty^ = e R ^,,' = ± ' J-00 fiw^ ' 



(A6) 



where (. . .)o denotes the average with respect to the equi- 
librium state determined by the Hamiltonian Tio, and Tk one obtains^ 
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where we have used the properties 



and 
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(A8) 



(A9) 



Since the electron-electron interaction ([8]) contains only forward scattering terms, all non-vanishing wire correlation 
functions must involve an even number of operators with a given chirality r. This yields r2 = —ri = ra, so that 



/^2^(x,t) = — 7+7_ 2^ A,; 



dtidt2dt3 ^ m'>l2i]3 ^ 
V,Vl,V2,V3=± r3 = ± 
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The term with ^3 = -I- can be shown to yield the same contribution as the term with 7^3 = —. To see this explicitly, 
one makes use of = — exploits Eqs. (|A8|) and (|A9|) . and renames variables according to —r], rji ~r]i 
(i = 1, 2, 3) and ti ^ t2. One can then write 
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where 



I7=± 

contains correlation functions of wire operators, while 



(A12) 



(A13) 



is a correlation function of the tip. These correlation functions are evaluated in App. [B] starting with Eq. (jB7p and 
(|B23|) . respectively. Inserting these results, one obtains 



SfjJ gje[(ys+Vt>-2yT)(ti-t2) + (Vs-VD)(ti-t3+t2-t3)]/2fig-2i[few+3^e(ys+V'D-2Vb)/?ii'w](xo-a;i) 



ril,V2,ri3=± 



(AM) 



X 



Kcl 



dx-'dtVw(x')(7^K[5,e(a;,i)a,,$(x',t')])r 



where, for any pair of bosonic operators A and B, the following definitions hold 

(A(tA)i3(tB))^^' - {{A{tA),B{tB)}) 
{AitA)B{tB)y'' = e{tA - tB){[A{tA), B{tB)]) 
{A{tA)B{tBT^^ = -e{tB - tA){[A{tA), B{tB)]) 



(A15) 
(A16) 
(A17) 



We now observe that the last term in Eq. (jA14p can be dropped. Indeed, since it depends neither on rfi nor on ti 
(i ~ 1,2,3), it can be singled out of the sums and integrals J dtf, the remaining sums and integrations yield a 
vanishing result, since the corresponding expression equals the term of order 7^ A of an expansion of 



Simple transformations of the integration variables of Eq. (jA14|) , and use of the relations 



Kcl 



dt {d^Q{x,t)'Prixo,0))o 



dt{d,Qix,mrixo,0))o 



4uw 



[1 + r sgn(a:: — xq) 



(A18) 

(A19) 
(A20) 



obtained from the correlation functions provided in App. [Cl yield 



1^2)^{x,t) 



IGTr^axa^ 



7+7 



- ^ A, / / dhdt 



X'^l QH{Vs+VD-2VT)(ti-t2) + iVs-VD){ti-t3+t2-t3)]/2h ^-2ilkw+g''e{Vs+VD-2VG)/hvvf]ixo-Xi) 



X J2 'n^2%F^"=''^(ti,t2)F:^^"=(ti-t2)b™,(tl,t2)b:[y^^,,^^(tl-t2) 

';i,I72,T)3 = ± 



(A21) 



X 



[Vi +V2- 2773 + sgn(a; - a;o)(77i - 772)] 
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Taking into account Eqs. (jBSp . (|B16p . (jB24[) . and (jB26|) . we now observe that upon reversal of Keldysh contour indices 
-Vi (« = 1,2,3), 



implying that in Eq. (jA2ip the contribution for 7/3 = — is conjugate to the one stemming from 773 = +. Thus 



(A22) 
(A23) 

(A24) 
(A25) 



7+7 



w 

CjJ gie[(Vs + VD-2VT)(tl-t2) + (Vs-VD)(tl+t2)]/2fig-2i[/cw+/e(Vs+VD-2VG)//i^'w](a;o-a;,) 



(A26) 



X E F^"^+(^i,*2)F^^^^(ti-i2)b;j^;'^=+,,(ii,t2)b^;:j=.(^^ 

77l,?72=± 

X [772 + 771 - 27717/2 + sgn(a; - a;o)(772 - 7]i)] 

The term sgn(x — xq) appearing in the last line is positive (negative) for a measurement point x located in the drain 
(source) lead. Recalling that the current can be written as in Eqs. p6|) . (|17p . it is easily seen that those terms 
that are multiplied by sgn(a; — a;o) yield /t/2, whereas the other ones yield Im- Inserting Eqs. (|B8p . (|B16p . (|B24p . 
and (|B26p into Eq. (|A26p . and changing to dimensionless integration variables = tv^/gL, the result (|69p is obtained. 



Similar procedures can be applied to evaluate the terms of order A^, A'^ and 7^. We find 

„3 



(A27) 



I^4x,t)^-^ XA.xJ dhdt.dt, ^i^2 3{w^r'|,,(ii,i2,t3)e-2^'^--(^'-^^)} 



(A28) 



and 



1^2 (2;, t) 



ev^VT 2 



r=± 



rii,ri2=± 



(A29) 



where 



^'iSj.ihM) = E (^K [j('')(:r,i)Z^w*V''^(-^-*i)*-^'^(^-^i)*-^'''^C^J-'^2)<''^^( 



(A30) 



W'j:,%{hMM) = E ^^(rK [j('')(x,0Ww*V''H^-*l)*-'^(^-*l)*-'''(^J-'^2)<''^^(x,,i2)p^ 



and 



^"^ZihM) = E (^K [/''n^,OWw*J^''^H^o,ii)^l''^n^o,t2) 



(A31) 
(A32) 



APPENDIX B: EVALUATION OF W AND 
T-FACTORS BY BOSONIZATION 



'Hq.w + T-Lf);T- For a non-interacting wire 7io,w = 'Wkin,Wj 



The Hamiltonian (jA3P of the interaction picture de- 
composes into commuting wire and tip parts, i.e., Tio = 
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and the wire correlation functions W introduced in 
Eqs. (Uni), (|M0)) . ((MT|1 and ((X32|) can be factorized 
into products of single-particle electron correlators using 
Wick's theorem. In this case the W's can be evaluated 
straightforwardly, and the results for the contributions 
(|XTT|) . (|M7l) . ([X28|) and (jMOl) to the current coincide 
with the corresponding terms of an expansion of the cur- 
rent obtained from the scattering matrix formalism. In 
the interacting case, however, ?io,w — ''^kin.w + 7if7i and 
Wick's theorem cannot be applied. In this appendix we 
evaluate the wire correlators W using the bosonization 
technique4i The wire field operators can be represented 
as 



(Bl) 



where the fields <I>± describe particle-hole excitations, 
and Klein factors represented as Majorana 

fermions.— Finally, aw is a cut-off length of order the 
lattice spacing. 

Introducing Eq. (|Bip into Eqs. ^ and (fT51) . one obtains 

Wo.w = ^ J^dx |:n2(x) + [d,<i>{x)f^ (B2) 

where $ = $+-(-$_ and 11 = —dx{^+ — are con- 
jugate bosonic fields, i.e., [^{x,t),Il{y,t)] = i5{x — y). 
Finally, 



1 

u 



for |a;| > L/2 

-1/2 (B3) 
for < L/2 ^ ' 

I 



is the inhomogeneous interaction parameter. Notice that 
< g < 1, where g = 1 describes the non- interacting 
case present in the leads. The limit g corresponds to 
strongly repulsive interaction. The wire current operator 
Eq. is expressed in terms of the dual field 6 = $+ — 
as 



I{x, t) = evw{dx&{x, t)) 
Further, with the help of the relation 

Pr[X,t) = ^ , 



(B4) 



(B5) 



the term ([6|) of the Hamiltonian can be written as 
1 /■+°° 

U^^ ^ —= dxiJ.w{x)dx^{x) . (B6) 

We start by discussing the derivation of 
W^3l','"(ti,t2,t3). Inserting Eqs. dH and 

(jBSP into Eq. ()Alip . one obtains 



W™^(tl,t2,t3 



(B7) 



1 



(27raw) 



2^ W 



Here 



Fr="^(ti - - h) - (Tk [4'^)(ti)^c^'^'fe)^^L"^''(t3)4'^^fe)] )^ (B8) 

e{t3 - ti)e{h - t2)f]n]2 + e{ti - h)e{t2 - h) - 6{t2 - h)e{h - ti)j]im - o{ti-ti)e{h - t2)mm 

accounts for the correlation function of fermionic Klein factors, whereas 

B^;'^?^,,(ti,t2,t3) = 



^^5J^^\x,t) 



Tk exp U -i ^ ^' / dxVw(x') J^"" ^ + j d^' 

\ I \ ,,' = ± ri"=± 



(B9) 



correlates bosonic vertex operators. Also, we have introduced the notation x = {x,t). The expression (|B9[) can 
straightforwardly be evaluated taking into account that for a functional 



C[J] = ^Tk jexp + ^ J dxj('')(x)i?('')(x)^ 1^ , 



where A and i? are linear combinations of bosonic operators, one haa^ 

5Ji^) (x) 





AS^'")(x) 


J=0 ^ 





(BIO) 



(Bll) 
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Furthermore, it can be shown that Ti.^^, i.e. the first term appearing in the exponent of Eq. (|B9p . simply yields a 
shift in the operators $± according to 



where the zero modes 



^ , 1 e[^s-VD+r(yg + FD)] ^ 
$o,r(a;,i) = — r t 



(B12) 



(B13) 



5'(Vs+Vb-2FG)^ fora;<-i/2 



40r ftww 



[5'(l^s + I^D - 2Fg) + r(Fs - Fd) ] : 



for \x\ < L/2 



{l + r)x-- 



L 



+ 9' {Vs + Vu-2Vg)- for x > L/2 



fulfill the equation 



$0,r(x) - $0,r(y) 



dx'/.w(x') [($,(x)a,<i>(x'))'-°' - ($,(y)9,$(x'))^ 



(B14) 



After lengthy but straightforward algebra one obtains 

B™,(tl,t2,t3) = 

= -2ie''''/''n^s(ti-t3)-Vb(t2-t3)-s^(Vs+VD-2VG)(a;o-a:i)/fw] ^';i';2»)3 U ~ U) 



-(0,e(x,t)$(.T„t3)>r - ^3(a.e(x,0$(a;„t2))r + ^ / dx'^iwix'){TK [dM^,t)d.'<^{^')])T } . (B15) 



where 



bl'^^^l,{h,h) = exp^-2^(TK 



(B16) 



exp{4^[7^w(eo;6:;n;T2) + iX^''^''-^(^o;e.;n;T2)]} . 



The correlation fimctions TZ^{^o;^i;Ti;T2) and 1^^'^^^ {£_o;$i;Ti;T2) are defined in App. \C\ [see Eqs. (|Cip - (jC6P ] 
and also given explicitly there in the zero temperature limit. The arguments = tiV^^/gL — tiU*^ and = Xj/L 
(j = 0, 1, 2) are dimensionless time and space variables. In deriving Eqs. (|B15|) and (jB16[) we have used the equalities 



and 



{A{tA)B{tB))'^-' = \ {A^''-'HtA)B'^''^\tB))=2^{A{tA)B{tB)) 

{A{tA)B{tB)y'' = I MA^''^\tA)B^''^\tB))^2\e{tA~tB)^{A{tA)B{tB)) 
r]A-VB=± 

{A{tA)B{tBV^ ^ I J2 VA{A^''-'\tA)B^'^^\tB)) = -2\6{tB -tA)^{A{tA)B{tB)) 



(B17) 

(B18) 
(B19) 
(B20) 



riA,riB=ii 

valid for any pair A and B of real Bose operators. 
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As far as the tip correlators T are concerned, see Eqs. (jA13|) . (|A29[) . and (jA27[) . Wick's theorem might be applied, 
since the tip is supposed to be non-interacting, and the use of bosonization is unnecessary. However, to have a uniform 
formalism and notation throughout the paper, we prefer to utilize a bosonized approach for the tip as well. The tip 
electron field and density are written as 



c{y) 



and 



c:[y)c{y) := , 



(B21) 



(B22) 



where ip{y) is a chiral (right-moving) boson field, and kt and ax are the Klein factor and cutoff length of the tip, 
respectively. By way of example, wc evaluate here the T- factor (jA13|l appearing in the calculation of I'^2x. Inserting 
Eqs. (|B2l1) and ((B2^ into Eqs. and ((IT3|) . one obtains 



^ 27raT 



F^^''=(il-t2)B^^^^^,(tl-i2) 



where, similar to the wire case. 



(B23) 



(B24) 



accounts for the correlation function of fermionic Klein factors, whereas the correlator of bosonic vertex operators 
reads 



\ ri' = ± ^ 



(B25) 



It is easily verified that the first term in the exponential function, which originates from the term (|12p in the Hamil- 
tonian, merely yields a time-dependent phase factor, so that 



where 



h^'^lliti-h) = exp -2^(rK 



93(''i)(0,ti)-¥>(''^)(0,i2) 



(B26) 



(B27) 



= exp{4^ [7^T(rl - t^) + iX^^^^Cn - t^)]} 
The functions 7^T(T) and J^^''" (r) are given in App. O see Eqs. ((CTT|) and ((CT2)) . 



APPENDIX C: CORRELATION FUNCTIONS 



This appendix collects properties of correlation functions appearing in Eqs. (|60|) and (j69p . as well as in Eqs. (jB16p 
and (jB27p . The transport properties of the wire are expressed in terms of the functions 

7^w(eo;6;rl;r2) = (Co; n) + (Co; 6; ^2) + 7^*+*- (Co; C.; n) + (^o; 6; ^2) (Ci) 

-7^*+*-(eo;eo;rl-r2)-7^*-*+(e.;6;0) 

ir'''-Meo;C.;ri;r2) = ^{M(ri)-r;i0(-ri)]X*^*'-(eo;C.;ri) + M(r2)-r720(-r2)]X*-*-(eo;C«;r2)} 



- [md{ri - T2) - md{r2 - Ti)] x*+*-(Co;eo;Ti - T2) 



(C2) 



where the functions 7?.*g*''(^; r) and X*'^*'' (^; r) are the real and imaginary parts, respectively, of the auto- 
correlation fimctions of the bosonic fields <I>r. Specifically 



<cf'(e;e';r) = 3?|^$,(x,0$,(x',o)-i[<i>2(x,i) + $2(x',o)]^ I 
^*"*'-(C;e';r) = 3{(<i>,(x,t)ci>,(x',o))o} 



(C3) 
(C4) 
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Likewise, the real and imaginary parts of the cross-correlation functions of fields with different chirality r read 



-R-"-"- (C; e-, r) - m.r{y, 0))o} 

I*"*-" (^; C'; r) - 3 o))o} 



(C5) 
(C6) 



Notice that the real part of the correlation functions of fields with the same chirality needs to be defined with an 
infrared regularization as in Eq. (jC3p . The above equations are given in terms of the dimensionless time and space 
variables r = tv^^/gL and ^ — x/L introduced previously. From the inhomogeneous Luttinger liquid model one 
obtains at zero temperature 



and 



62n ^ — ' a 



[g + g +2r-) 2^ p' 'In ^ 



+ ig-g-') E In 



™SZodd 



+ ln 



21 2 



[a2^ + (2^ + m)2][a2^ + (2C' + m) 



I JTieZcvon 



+ (5 + 5 ' + 2r) E 



I m 



arctan 



r + + m 

aw 



(C7) 



(C8) 



and 



+ {g-g 



arctan | ''" ^ '^"^ | + arctan ' ^ 



aw 



aw 



= -^((5-5-) E 



?7ieZevcn 



1^ aw + (^ + ^r+m)^ ^ ^ 1^ a^ + (r - - m)^ 



"w + 



+ (3 + 3 + 2r-) > p' ' In = — — — — • — - 

[aw + (^fl + m 



"iGZodd 

+ (3 + 3"') E ''""'in 



,21 2 



[a2^ + (2e + m)2][a2v + (2C'+m)2] 



4 E 



|m I 



in 



+ In 



\a% + (2C + 1 + mf\ [a^ + (2g' + 1 + m)^] 
("w + "^^)^ 

[aj + (2g - 1 + mf\ \a%, + (2^' - 1 + m)^' 
(a2^ + m2)2 



(C9) 



167r 



-5-^) E 



?7ieZcvcn 



arctan 



T + 5r + 



arctan 



(5 + 5 ^ - 2r-) E P''"' arctan 



aw 

r + + m 
aw 



aw 



(CIO) 



+ (9 + ^ + 27^') p'™' arctan 

"iSZodd 



Here we have introduced = C^S'i C-R = C + S'; and the dimensionless cutoff length aw = o-w/gL, as well as the Andreev-type 
reflection coeflicient p = (1 — 5)/(l + 5)- 
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The correlation functions for the non-interacting tip can directly be obtained from the above results. The tip is described 
by a single chiral mode, and we need the correlation function only for coordinates at the injection point y = 0. From Eqs. (|C7[) 
and (|C8|) we find for ^ = 5' = by taking the limit <? — > 1 and replacing aw by ot 

7^T(r) = 7^rg(0;0;r) = 7l*|*+(0;0;r)| ^_i.ln^l_^ (Cll) 

J^i^^^M = F:[,i''=^(r)X*'*'(0;0;r) = F:[,i''2(r)2:*+*+(0;0;r)| = _ r;26>(T) arctan f-^\ . (C12) 
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